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We use a functional integral formalism developed earlier for the pure Luttinger liquid (LL) to 
find an exact representation for the electron Green function of the LL in the presence of a single 
backscattering impurity. This allows us to reproduce results (well known from the bosonization 
techniques) for the suppression of the electron local density of states (LDoS) at the position of the 
impurity and for the Friedel oscillations at finite temperature. In addition, we have extracted from 
the exact representation an analytic dependence of LDoS on the distance from the impurity and 
shown how it crosses over to that for the pure LL. 



I. INTRODUCTION 

One of the first exactly soluble models in the prob- 
lem of strongly correlated electrons was formulated in 
one dimension in the seminal papers of Tomonaga^ and 
Luttinger— and solved by Mattis and Lieb'^. Consider- 
able further contributions to understanding of generic 
properties of the ID electron liquid have been made in 
paper a^i^i^ . In particular, Haldane^i has coined the notion 
of a Luttinger liquid (LL), stressing the generic proper- 
ties of the Luttinger model for ID Fermi systems, and 
has formulated fundamentals of a modern bosonization 
technique as one of the most elegant ways for solving the 
problem. In this approach, the Fermi creation and anni- 
hilation operators are explicitly represented in terms of 
Bose operators and a 4-fermionic Hamiltonian is eventu- 
ally diagonalized in the bosonic representation. 

There exists an alternative way to bosonize the prob- 
lem (sometimes called 'functional bosonization') which 
was suggested ini and found its further elaboration 
in papersSiSiifl. In this paper, we will apply a func- 
tional method developed earlier by one of usAS for 
the treatment of a single-impurity problem in the Lut- 
tinger model. The essence of the method is in the de- 
coupling of the four-fermion interaction by the stan- 
dard Hubbard-Stratonovich transformation, typical for 
higher-dimensional problems, and a subsequent elimina- 
tion of a mixed fermion-boson term in the action by a 
gauge transformation which is exact for the pure ID Lut- 
tinger model and gives a convenient starting point for 
including a single backscattering impurity. 

The problem of the impurity in the LL is, of course, not 
new. However, it was a subject of some controversy. Al- 
though it is generally acceptedii*i2iH that the local den- 
sity of states (LDoS) at the impurity site is suppressed in 
agreement with a generic consideration due to Kane and 
Fisher fi4 an alternative prediction of its enhancement as 
also been made Although our method is rather similar 
to that of Ref. 15, we show t hat its consis tent application 
leads to the results of Refs. Illll^l.-lld Moreover, we 
obtain a general expression for the Green function appli- 
cable in a wide range of temporal and spatial arguments 
that contains not only the suppression of the LDoS at 
the impurity cite, but also describes its dependence on 



the distance from the impurity, as well as the Friedel 
oscillationsSiiSiii. Such a generality show the suggested 
formalism to be quite powerful and useful. 



II. EXACT FUNCTIONAL REPRESENTATION 
FOR GREEN FUNCTION 

The Hamiltonian of Luttinger liquid with one 
backscattering impurity can be written as: 

d_ 

' dx ' 



i / dxdx'i:l{x)i^l,{x')Voix - xyj^,{xyj^{x) 



+ v^Jdx\{x) ijl_{x)i;-{x) +4!l{x)tp+{x) (1) 

Here V'± ip± are the standard creation and anni- 
hilation operators for left- and right-moving electrons 
(-0 = ip+e^Pp^ + ip-e~'^Pi'^), T] ~ ± and the summation 
over repeated indices is implied; Vq is a bare electron- 
electron interaction. The impurity potential is given by 
WpA(x) — VpXu{x), where u{x) is some form-factor of the 
impurity and A <g; 1 is its strength. 

The temperature Green functions of the Hamiltonian 
can be represented by the functional integral: 



where e^(.,.),Z.^/e-[^lp^, and the action 



(2) 



(3) 



Here the integration over the imaginary time r goes from 
to P = 1/T, and the "classical" Hamiltonian H is ob- 
tained from H, Eq. (^), by substituting the Grassmann 

fields ^^'^ V'(f) for ip'^ and 

The Green function of real electrons is given as a sum 
of all Gjj^' taken with appropriate phase factors: 



(4) 
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The functional bosonization is achieved via the stan- 
dard Hubbard-Stratonovich transformation decouphng 
the 4-fermion term in the action. Introducing a new clas- 
sical (bosonic) field (p we arrive at the action 



-V, I d^x{x) [r+io^^-io + r-ioi^+io] (5) 

Here Vq^^ is a function inverse to Vq in the operator sense, 
whose T dependence is just (5(t— r'), and the 'chiral deriva- 
tives' dri are defined by 

d+ = 2dz = 9r - iv^dx , d- = 2dz = 9t + iv^dx 

In order to eliminate the mixed term in the chiral 
derivative in Eq. we apply the gauge transformation 

V-.^V.e^'", r^^r^e-^'"^, d^e, = ^, (6) 

where 9^, is a complex function which depends on the 
field (/>. Since ((> is real, 



This transformation produces^*^ the Jacobian J, derived 
explicitly in the Appendix, 



(7) 



where 11 is the polarization operator whose Fourier trans- 
form is given by 



n(g,r!) 



1 



9 ? 



(8) 



where 17 — 2TrnT is a bosonic Matsubara frequency, 

Therefore, after the transformation © the interaction 
in the action ||SJ) becomes RPA-screened as expected, 
Vq V^^ = Vq'^ + n, and can be split into the sum 

S[i:, 0] = Sfii:] + S4(j)] + Si^pii:, 0]: 

s^ = ^ fd^demv-'i^-am 



(9) 

The Green function Q can be represented as the func- 
tional average over the fermionic and bosonic fields with 
the weight Sh + Sf. 



(10) 



The bosonic field (j) enters the pre-exponential factor in 
Eq. pUI) only implicitly, via 9{4)). Before proceeding fur- 
ther it is convenient to work out correlation functions 
of 9 which follow straightforwardly from {4'{£,)4>{0) ^ — 
V{^-C) and Eq. ©: 



, , , , 1 I sinfz — z')\^ 



2 |sin(Zp-2:;) 



(11) 



2 sin(zp-z') 



where 

= 7rr(T + ix/v^) 

Viq = 0) 



V = { 1 



TTWp 



z = ttT{t + ix/v) 

V, (12) 
5 = — ■ 



Here we assumed that the Fourier transform of the 
forward-scattering pair interaction only weakly depends 
on momentum, i.e. y(g^2pp) « V{q — 0). 

Now we reduce the partition function Z\ = ^^e~'^'"p)) 
in Eq. to the product of fermionic and bosonic in- 
tegrals. This can be done for an arbitrary scattering 
potential A(a;). On expanding e~'^''"p and keeping only 
the terms with equal numbers of and ip-i- (as well as 
of -01 and we obtain 



oo 2n ^ /• 

= E 7^ n Kfed^I.A(:rfe)A(4) 

71 

X (11 V';(a)^^(a)C(a)^+(ei.))^ (13) 

n 

x(exp [2^(02(a)-^2(efc))])^; 



k=l 



Carrying out the bosonic average with the help of formu- 
lae (|11() . we find 

(exp [2Y.m^) - 02{m), - J!.rL ^ (14) 



k=l 



where 

a ~ T/e^ < 1 
is a cutoff parameter, and Pn{z) is given by 



n sin(2j - 2j)sin(z,; - zj) 

Pn{^) - — ^ ; (15) 



Y\ sin(^ 



The parameters z and Zp , entering with the appropriate 
indices Eq. itHll . are defined by Eq. ifT^ . 
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The fermionic average in Eq. (|13|l is independent for 
the left- and right-moving electrons and yields 

" / T \" 1 



fc=l 



Applying the Cauchy formulaic, 



det 



sin(zi - z'A 



one sees that the fermionic average cancels |P„(zp)p 
in the denominator of Eq. (|14|l so that 



2n 



71 „ 

X n K/cda A(a;fc)A*(4) |P„(z)P^ . (17) 
k=i'' 

The result above has been obtained by formally calculat- 
ing both fermionic and bosonic integrals for the partition 
function in Eq. H10|l . Now we 're-bosonize' this expres- 
sion by presenting it as a result of the integration over a 
new bosonic field 6: 



-SA[e]\ 
/ 



(18) 



where 



T 



Sxm = /deA(a;)cose(C) (19) 

a J 

and (...)o average is defined with the action Sq: 

Som = ^ / [i^rOf + v^d^ef] (20) 

To verify the validity of the representation p8|l -H19 |l one 
needs to expand the exponent in Eq. H18|) using the fact 
that the pair correlation function of O with the action Sq 
is given (with a proper regularization) by 

Goita ^ (e(O0(C'))o = -25ln| sin(z - z')\ . (21) 

The sum resulting from such an expansion coincides with 
that in Eq. H17() . We remind that here ^ stands as a 
short hand for x, r with the appropriate indices, while 
z = TTT{T + ix/v), Eq. (P|l . 

On repeating the steps outlined in Eqs. (|13|I - H17() . we 
obtain the following 0-representation for the Green func- 
tion of Eg ifTHI: 



T 



^e,,,(e,c'), (22) 



2v^a^-^/^a ^s,m{z-z')\'^ 
where we introduced an auxiliary function ^jjjj' (Cj C')- 
Q,A^^O = Z-,' (c^««)-^«(«')e-^^[«-^l)^ . (23) 



The COS0 term in action Eq. H19|l . is now shifted by 
the phase factor x(^i) 

Sx[Q-x] = -- /dCiA(xi)cos{e(6)-x(6)} (24) 

where x('?i) parametrically depends on the arguments of 
the Green function 122II . 

X(xi,Ti) = x(2^i) = arg— -. (25) 

sm(zi — z') 

Finally, the sign factor s^^/ in Eq. (|22|l is defined (with 
■q = ±1) by 

Sijjj' {z — z') = exp [\i{ri -\- ij') arg sin(z - z')] . 

In the particular = 0, the representa- 

tion (|22|l -123 () coincides with that obtained by Oreg and 
Finkcl'stein^^, although we will show below that it leads 
to the well known resullii^ for the density of states rather 
than that declared in RefE 



III. SELF-CONSISTENT HARMONIC 
APPROXIMATION 

The representation H22() - (|23|l is exact and further cal- 
culations are only possible after some approximations. 
First, we assume the impurity to be point-like, X{x) = 
XS{x), and weak, i.e. A ^ 1. This ensures the spectrum 
linearization to be valid in the presence of the impurity 
and all the relevant energy scales to be small compared 
to the Fermi energy. 

It is well known that even such a "weak" impurity 
leads to strong changes to the single-electron density of 
states in its vicinity. It's influence is perturbative only 
in the high-temperature limit, A <C a oc T/ef- In the 
low-temperature regime, A 3> a, a non-perturbative ap- 
proach is required. In the present context, it can be 
developed within the so called self-consistent harmonic 
approximation (SCHA)^^. It is based on the fact that 
for a ^ A ^ 1, the deviation of 8 from x('''i) in the 
action (|23|l is prohibitive so that their difference can be 
presented as small quadratic fluctuations around one of 
the minima of cos(0 — x)- By minimizing the difference 
between actual cosine-shaped potential and its quadratic 
flt (the Feynmann- Vernon variational principle) , one sub- 
stitutes the exact cosine potential by the harmonic one 
thus reducing the action H24|l to the following one: 



AT f 

Sk[Q -x] = ^JdT, [6(0, n) - x(0,Ti)]' , 



(26) 



where the new 'impurity strength' A is chosen to pro- 
vide the best flt to the real potential. As a result, 
one obtains^^ the renormalized 'self-consistent' impurity 
strength A as A = A^^ . 

Now the action is quadratic in 0, 



AT 
2^ 



J dre2(0,r) 



(27) 
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so the integral (|23() for Qjjjj' {£,, CO is reduced to calculating 
the averages with the action H27|l of linear in Q terms in 
the exponent, using the standard formulae of the type 



where (. . .)q stand for the functional averaging with the 
action 12 7|) . The integration thus yields 



It should be stressed that both S and $ depend on the 
'observation points' C and C' via the appropriate depen- 
dence of the parameter x, Eq- ij^ - AH these functions 
can be calculated for arbitrary ^ and C'- However, since 
we are only interested in the local density of states (at 
an arbitrary distance from the impurity) and Friedel os- 
cillations, it is sufhcient to consider x — x' case only; r' 
for convenience is set to zero so that from now on we use 
C — {x. t) and C' = (a;, 0). Then we obtain 



lnew(C,C') =^G(e,e) + ^G(e',e') - ^G(C,C') 



(28) 



All the terms above can be expressed via the pair corre- 
lation function of the auxiliary bosons Q defined by 



G(C,O-(0(Oe(n> 



e ■ 



(29) 



In the absence of the A term in action (|27|l . G reduces 
to the standard bosonic Green function Go, Eq. (|21|l . 
The full Green function G in the presence of impurity is 
straightforward to find in the Matsubara frequency rep- 
resentation, 

G(C;e') = r^G(x,x';c.)e-'"(^-^'), 

where it is expressed via Gq{x^x' ]uj) as follows: 

AT Go{x,Q-uo)G^{Q,x'-uo) 



G{x, x' — Gq{x, x'; lo) 



27r.g 



- \x — x' I 



a 1 + M:Go(0,0;w) 

g-M(|a;| + b'|) 



27rgAT 



1 



(30) 



Exponentiating the denominator of the second term 
above by 1/D — 2 Jdse^^'^'*, we obtain G in the x,t 
representation: 



G(e,e')-Go(e,e') 



(31) 



J dse' 





2s 



In 



ttT t — t'- 



\x\ + \x'\\ as 
V J gA 



The impurity-induced terms in Eq. (|28|l . $ and S, 
which are respectively linear and quadratic in the fac- 
tor X, Eq. (|25(l . result from the averaging of the first and 
zeroth order terms in O arising from Eq. I|2t)|) . 



AT 
2^ 



dTldT2 



rAT 



G(0,Ti;0,T2) 



S{ti - Ta) x(n)x(T2) 
AT 



(32) 



2a 



dTi[aG(C;0,ri)-a'G(e';0,ri)]x(Ti) 

(33) 



2a 



AT 

1 G(0,0;w) 

a 



^x(0,c.)[aG(a;,0;c.)e-'"^- 



AT^ 
' ~2^ 

a'G(0,x;w)] 



(34) 



(35) 



Substituting the Fourier transform of Eq. H25|l . 



, , 171 I .1 1^1 

X{0,(jJ) = — sgnxe ' ' " 



(36) 



and Eq. into Eqs. and (123, we carry out the 
Matsubara summation to obtain S and $ in the same 
representation as follows: 



C30 

S(x;T)-^ydse-2^ In 



sin ttTt 



smh'[fi + 2nT 



\x\ 



(37) 



7] — J] 

$^,,'(a;;r) = ——sgnx 



3m In sin ttT [ t + 2i — 

2 \ V 



Now we have all the ingredients to find the Green func- 
tion Using Eq. |(2SJl, we express Q{£) in terms of G, 
$, and S as follows: 



T 



Q,2g 



I sin ttTtI 2s 



.f.-^(.OQ-\G(x,x-r=0) 



X sgnr e 



cos[2ppX + $+_(e)] (38) 



The appropriate asymptotic behavior of G, $, and S is 
determined from formulae H31|l . (|37|l . In region 5|a::|Pp <^ 
, ItI^f <C A""'^ the impurity can be described pertur- 
batively and its presence results only in a small correction 
to the Green function of the pure LL. So we restrict our 
attention to G{£) = Q{x, x; r) outside of the perturbative 
region on the {x, r) plane. Introducing convenient dimen- 
sionless notation x — gp^\x\, f — Sf\t\ and keeping only 
the leading contributions, we find 
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' 2glnA-\ i<A-i;r = 



G{x, x;t) = < 



2.9 In ^ ^, a;>A"\r 

sm ar 



0, 



1 , sinaf . , 

In A^i < a; < r 



2g ax 
1 

23 



-fx t) = < 1 , Asinaf , , 



0, 



X » A^^f 



We do not write asymptotics for $_| (x^t) explicitly since the phase does not make any real impact. The above 

expressions enable us to find the Green function, Eq. H38() . for any distance x from the impurity: 



g{x, x; r) = < 



(sin af) ! 



T a 



(sin af) 9 



. ''^F (sinaf) 2 ( 



sgnr — cos(2ppa; + <i>^ ) 

sgnr — cos(2ppX + $_| ) 

/ sinaf \^ ,„ 
■gnr - -— — r cos(2p a; + 9 
\smh ax J 



X < A^^ < f 



A^i < i < f 



r 



To begin with, the general expression for the Green 
function, Eq. H38|) . allows us to extract a well known re- 
sult for the Friedel oscillations&iSiii 

T 1 

Ap{x > (ApJ-i) = J— — -— gCos2ppX; 

^ I V \ 

(39) 

For X in the range (ppA)~^ <C x ^ (^t being the 
thermal dephasing length) the amplitude of the Friedel 
oscillations decreases as ~ |a;|~^, and for larger x it is 
exponentially suppressed at distances exceeding ix/'ig. 

More interesting is the local density of states (LDoS) 
which is defined via an appropriate analytical continua- 
tion of the Fourier transform of Q: 



v{x^e) = 3m / d7 

71- J 



Using asymptotics for the Green function, we get: 



z^(x, s) 



en 



lA-Ml-^), 



£21,9 



X < A^i < e^^ 
A^i < X < e"^ 



X > A"\£"i 



(40) 



The three regions above with different behavior of LDoS 
are sketched in Fig.l. The first line in H4Q(I describes 
LDoS in the vicinity of impurity, in full correspondence 
with the original results of Kane and Fisher^'* obtained 



for the LDoS at x = 0, i.e. exactly at the impurity. In 
addition, we specify here a dependence of the LDoS on 
the strength of the impurity (we want to remind that A = 

1 _ l + g 

A 1-9 , which leads to LDoS being proportional to A ) 
The region of applicability of this result corresponds to 
the diagonally hatched section in Fig. 1. 

The analytic dependence of the LDoS on the distance 
from the impurity given by the second line in Eq. H40f) de- 



A 




hi ii,K 



r 3 is c \ er 



HE 



FIG. 1: regions with different behavior of local density of 
states i'{x,e); x = gpp\x\, f = £f|t| 
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scribes a new result for the crossover from the impurity- 
induced dip in the LDoS to the bulk behavior. So, as 
a function of x at a fixed e, ^{x) remains constant for 
X < A~^, then it increases until x becomes of order e^^ 
(which happens before x overtakes the thermal dephasing 
length, £t — Vp/nT), where the LDoS reaches its bulk 
(x-independent) value, given by the last last line in (I40|l . 
This corresponds to a well-known result (see Ref. 20 for 
reviews) for the LDoS in the homogenous Luttinger liq- 
uid. Its region of applicability is horizontally hatched in 
Fig. 1 (although Eq. I|40|l does not describe the pertur- 
bative region i,e^^ ^ the presence of the impu- 

rity results there only in small corrections to bulk be- 
havior). The crossover between the impurity-dominated 
and bulk values is governed by a universal power law: 

In conclusion, we have demonstrated how powerful is 
the method of "functional bosonization" developed here 
for the description of a single back-scattering impurity in 
the Luttinger liquid. It allows one to get within a single 
formalism the results both for the local density of states 
at an arbitrary distance from the impurity and for the 
Friedel oscillations. 



r„ with n external lines corresponding to <?!)'s, each loop 
being built of the n Green functions g± , respectively: 



/n n 
Y[dxkdTk r(f)(zp^; ...;Zp„) ]^0(a;^,T,), 

(A2) 



fc=i 



p(a) 



with the boundary condition Zp^_^_^ = Zp^. Substituting 
(7^ from Eq. l|Aip . one finds 



r+fz ■ 



Si+1 



One can rewrite the symmetric part of this vertex, which 
contributes into the integral in Eq. ljA2p . as follows: 



n 

k=l 



Sk ■ 



Acknowledgments 

We gratefully acknowledge support by the EPSRC 
grants GR/R95432 and GR/S29386. 



APPENDIX A: JACOBIAN 

The Jacobian of the gauge transformation Eq.® can 
be written as 



In J[0] = Trln 



d„ 



00 ^ 



7?— ± n—1 



where the Green function of non-interacting right- or left- 
moving electrons, obeying — 1, are given by 



1 



T 



2v^ sin(2p - 2^) 



(Al) 



where Zp is given by Eq. H12|l . The n-th order term in 
(j) is the sum of two vertices made of the loops and 



where An is an absolutely anti-symmetric polynomial 
built on n variables Si. By power counting, its order 
should be n{n — 3)/2. On the other hand, the minimal 
possible order of an absolutely anti-symmetric polyno- 
mial of n variables is n{n + l)/2, as follows from the 
fact that the powers of different variables should be dif- 
ferent for any monomial in order the entire polynomial 
to be anti-symmetric. The two inequalities can only be 
satisfied for n < 2 so that An>2 = 0. Therefore, all 
loops containing more than two external lines are zero. 
(Such an observation was first made by Dzyaloshinskii 
and Larkini within diagrammatic techniques). 

Therefore, we are left with the contributions from the 
loops with one or two external lines. The loop with one 
external line is proportional to the zero-momentum mode 
of the Coulomb interaction and is cancelled, as always, 
due to electroneutrality. The loop with two external lines 
is just the standard polarization operator in the random- 
phase approximation (exact for the LL), given in {q,fl)- 
representation by Eq. © so that the Jacobian is reduced 
to that in Eq. ((TJ in the main text. 
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